Charge mobility of discotic mesophases: a multiscale quantum/classical study 
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A correlation is established between the molecular structure and charge mobility of discotic 
mesophases of hexabenzocoronene derivatives by combining electronic structure calculations, Molec- 
ular Dynamics, and kinetic Monte Carlo simulations. It is demonstrated that this multiscale ap- 
proach can provide an accurate ab-initio description of charge transport in organic materials. 



Organic photovoltaic (OPV) and other organic opto- 
electronic devices rely upon blend films of two material 
phases, combining the attributes of a large interfacial 
area for charge separation and recombination and effi- 
cient vertical transport paths. Transport in organic elec- 
tronic materials proceeds by incoherent hopping and de- 
pends upon local molecular ordering as well as on the 
presence of percolation paths, and therefore is an intrin- 
sically multiscale property. Improved mobilities are cru- 
cial to improved OPV device performance, and in turn 
require control of structure on both the intermolecular 
and macroscopic length scales. 

Discotic liquid crystals could be ideal materials for 
OPV as they offer optimal design possibilities for both 
the charge mobility and the necessary underlying meso- 
scopic morphology. Discotic thermotropic liquid crystals 
are formed by flat molecules with a central aromatic core 
and aliphatic side chains [lj. They can form columnar 
phases, where the molecules stack on top of each other 
into columns, which then arrange in a regular lattice. 
Along the stacks of aromatic cores in the column one 
observes one dimensional charge transport 0, Per- 
pendicular to the column axis the charges have to tunnel 
through the insulating side chains, resulting in a much 
reduced mobility. Therefore the discotics in a columnar 
phase act as nano-wires and are a promising alternative 
to conjugated polymers, where charge mobility in devices 
is often limited by inefficient inter-chain hops. Columnar 
mesophases could, additionally, provide for laterally seg- 
regated phases of donor and acceptor molecules with high 
interface areas and efficient percolation pathways. The 
reported high charge mobilities [2j would guarantee ef- 
ficient operation. However, the spatial arrangement of 
stacks is never perfect: the columns can be misaligned, 
tilted, or form various types of topological defects, which 
are metastable or even stable at ambient conditions. In 
a recent molecular dynamics (MD) study we have shown 
that the degree of order in discotic stacks of hexaben- 
zocoronene (HBC) derivatives is sensitive to the type 
of attached side chains 0]. The local alignment of the 
molecules in columns and the global arrangement of the 
columns in the mesophase, are thus both sensitive to 
molecular architecture. 



A key question is the influence of local ordering in dis- 
cotic mesophases on charge transport. Computer sim- 
ulation is a powerful tool for probing the influence of 
molecular structure on charge transport. However, the 
task is challenging because different theoretical methods 
are needed to describe different length- and time-scale 
phenomena. Atomistic simulations are needed for local 
molecular arrangements quantum chemical calcula- 
tions arc needed for the electron transfer mechanisms 
(and interaction with the electrodes) [Bj], and stochas- 
tic or rate-equation dynamic methods for simulation of 
charge dynamics. In this paper, we present the first sta- 
tistical mechanics evaluation of the charge mobility based 
on quantum chemical and atomistic MD methods. We 
employ such a three level approach to explain the side 
chain dependence of charge mobility for differently sub- 
stituted HBC derivatives. 

Previous studies of liquid crystalline triphenylene 
derivatives Q have established that polaron hopping is 
an appropriate description of microscopic charge trans- 
port in these materials. To use charge transport theory as 
a predictive, rather than descriptive tool ab initio calcu- 
lations of the rates for charge hopping between molecules 
are needed. To this end we use the Marcus-Hush formal- 
ism |7| for the charge transfer rate u>ij 
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where Jjy is the transfer integral for electron or hole 
transfer, AGij the difference in free energy between the 
initial and final states, A the reorganization energy, h 
Planck's constant, k Boltzmann's constant and T the 
temperature. Jy and A can be computed using quantum 
chemical methods jB[. Note that Jy is highly sensitive 
to the relative position and orientation of the molecules 
involved, which will be determined, in turn, using Molec- 
ular Dynamics. If entropic effects and energetic disorder 
are ignored, AGy = F ■ dij, where d^ is the vector be- 
tween the center of the molecules involved in charge hop- 
ping and F the electric field. In this study we choose to 
neglect site energy differences due to electrostatic inter- 
actions between the molecules, because such differences 
vanish when the conjugated cores are parallel. Further- 
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FIG. 1: Stick diagram of HBC derivatives with different side 
chains. We studied alkyl chains of different lengths, C n , 
with n = 10, 12, 14, 16; branched side chains, Cio-e; and the 
dodecylphenyl-substituted side chains PhCi2. 



more, the exclusion of site energy disorder allows us to 
focus exclusively on the effect of configurational disorder 
on charge transport, through disorder induced variations 
in Jij ■ 

Once all u>ij are known, charge dynamics can be simu- 
lated by Monte Carlo (MC) methods [H or by a Master 
Equation (ME) approach [9(. Though several aspects of 
the problem are treated in the literature [nj , a com- 
prehensive study on large columnar systems does not yet 
exist. In organic molecular crystals, Deng and cowork- 
ers [HI showed that mobilities in pentacene can be repro- 
duced by mixed quantum chemical and molecular dynam- 
ics studies. Here we use a similar combination of methods 
to achieve a truly ab-initio description of charge mobil- 
ity in HBC: (i) molecular dynamics simulations (MD) 
are performed on discotic mesophascs of different HBC 
derivatives, (ii) positions and orientations of molecules 
are used to calculate the transfer integral J for all neigh- 
boring pairs of molecules in a column, and (iii) charge 
dynamics are determined using both Master Equation 
and kinetic MC methods. 

All studied HBC derivatives are shown in Fig. [T] The 
molecules consist of a flat aromatic core and six side 
chains. Experimental data on their structural proper- 
ties 0, EH, 0, EH \iBL as well as microwave conductivity 
(PR-TRMC) data [2, E3 are available. In some cases 
time-of-flight (ToF) mobility data are also available [l8l |. 

Molecular dynamics. In our simulation we adopt the 
united atom approach for the side chains, and consider 
explicitly only the hydrogen atoms belonging to the aro- 
matic rings of the central core. Details of the force 
field and the MD setup can be found in Ref. 4. The 
molecules were arranged in 16 columns with 10, 20, 60, 
or 100 molecules in a column. Production runs of 100 ns 
were performed at constant pressure of 0.1 MPa and tem- 
perature T — 300 K fixed using the Berendsen method 
with anisotropic pressure coupling. System configura- 
tions were saved every 10 4 MD steps (the timestep was 
2 fs) , yielding 200 snapshots for each run. 

Quantum Chemical Calculations. For each MD snap- 



shot, the overlap integrals Jij between intracolumnar 
nearest neighbors i and j were calculated. To calculate 
Jij , the aromatic core of each molecule, as output by the 
MD simulation, is replaced by a rigid copy of the energy 
minimized configuration, with the same axial and tor- 
sional orientation as the MD result. Because of molec- 
ular symmetry, the highest occupied molecular orbital 
(HOMO) and the lowest unoccupied molecular orbital 
(LUMO) of HBC are doubly occupied. Thus we have 
to calculate four transfer integrals for hole and four for 
electron transport [Io| . For each component we used an 
adaptation of the ZINDO method, because of its speed 
and its insensitivity to polarization effects 19]. The effec- 
tive J values are then taken as the root-mean-square of 
the four HOMO transfer integrals for holes and of the four 
LUMO transfer integrals for electrons [2(J. The inner- 
sphere reorganization energy for cation and anion radi- 
cals was calculated using unrestricted wavefunctions, the 
B3LYP functional and 6-31g** basis set. We found reor- 
ganization energies A of 0.13 eV for cations and 0.11 eV 
for anions. The outer-sphere contribution was neglected. 

Kinetic MC and Master Equation (ME) simulation. 
For kinetic MC, columns of molecules produced by MD 
are stacked periodically to produce a stack 1 /im thick. A 
uniform field is applied along the stack, a charge carrier 
representing either an electron or a hole is introduced 
near one end, and the charge drift is simulated by a 
continuous-time random walk algorithm. The method 
used for MC simulation was similar to that in Ref. [2l[ 
with adjustments for a disordered lattice. The charge 
mobility /i is obtained from the transit time t tr of the 
simulated transient via fi = L/ Ftt r where L is the stack 
thickness, ftr is taken as the point of intersection of the 
two asymptotes to the simulated photocurrent plotted on 
a log-log plot, as would be done in a ToF experiment. 

To compare these to steady-state mobilities along the 
column, we solved the linearized Master Equation [9( for 
continuous boundary conditions. The probability Pj of a 
charge being on site i obeys the rate equation: 

dPi/dt = [uijPj - WjiPi] , (2) 

where the sum is over neighbors j =fi i. fi is 
then obtained from the charge velocity v = fiF — 
Ylijtj i^jiPi ~ UijPj) {rj—Ti), where r, is the coordinate 
of site i in the direction of F. 

Representative snapshots of two of the systems are 
shown in Fig. [2J For C12, the aromatic cores are well 
aligned and the columns are evenly spaced whilst for 
C10-6 the branched side chains lead to a more disor- 
dered columnar structure, as studied in detail in Ref. 4. 
Here we mention only the parameters relevant for charge 
transport, namely, the intracolumnar molecular separa- 
tion h and nematic order parameter S. The ordering for 
C10 - Ci6 side chains is almost perfect, while PI1C12 and 
C10-6 are slightly more disordered. For the bulkier C10-6 
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FIG. 2: MD simulation snapshots of columns of 10 HBC 
molecules with (a) C12 and (b) C10-6 side chains at 300 K. 
Columns are pre-arranged on a rectangular lattice. 



we find numerous defects in the columnar arrangement. 
The intracolumnar separations are practically the same 
for all derivatives, as given in table HI 

To study the effect of different side chains on mobility, 
we performed ToF simulations of photocurrent transients 
on columns constructed from copies of each MD snapshot 
and then averaged over results for different snapshots. 
This sampling method is analogous to the experimental 
ToF situation where charges are photogenerated in paral- 
lel columns of different configuration, and the measured 
transient is the sum of displacement currents from all 
columns. For the same F, t tr varies significantly from 
snapshot to snapshot, i.e. the mobility is sensitive to the 
particular orientational and positional configurations of 
molecules in the column. Therefore, detailed tests were 
performed to determine the number of different columnar 
configurations as well as the length of the stacks needed 
to obtain reproducible results. Tests showed that 100 dif- 
ferent configurations are sufficient, and for 10 molecules 
per column the finite size error is less than 10%. 

The results of ToF simulations are shown in Fig. [3^,. 
Comparison of the shape of the simulated transients for 
C10-C16 and C10-6 a t the same F shows that the disor- 
dered C10-6 structures lead to relatively dispersive tran- 
sients while those for C10-C16 lead to the non-dispersive 
transients typical of ordered materials [2^ |. In our ap- 
proach, where energetic disorder is neglected, disorder 
in charge transport can only arise from disorder in the 
nearest-neighbour transfer integrals J, which in turn 
arises from disorder in intermolecular separation and ori- 
entation. Figure[3jD compares distributions of log I J| 2 . As 
the data clearly show, the disordered compound leads to 
a much wider distribution of J leading to a wider distri- 
bution of rates. In particular, the low J peak represents 
bottlenecks at discontinuities in the stacks, which act as 
traps in the one-dimensional transport. 

Finally, we compare our simulations to PR-TRMC 
measurements. PR-TRMC probes the sum of high fre- 
quency photoinduced conductivities due to both types of 
charge carriers, and therefore probes the fastest contri- 
butions to charge transport, sensitive only to the local 
disorder within well ordered domains. PR-TRMC mobil- 
ities are therefore appropriate for comparison with our 
simulation since the column sizes we are using are too 
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0.2 [23] 
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TABLE I: Electron (e) and hole (h) mobilities (crr^V^s -1 ) 
of different compounds calculated using time-of-fiight (ToF) 
and master equation (ME) methods, in comparison with ex- 
perimentally measured PR-TRMC mobilities. Also shown are 
the nematic order parameter S and the average vertical sep- 
aration between cores of HBC molecules h (nm). 



small to include topological defects. ToF mobility data 
would, in contrast, be heavily dependent on defects and 
on fluctuations with wavelengths longer than the column 
size we have used. The simulated mobilities are compared 
with experimental data in table [J and in figure [5J;. The 
systematic overestimate of /x by the ToF compared to the 
ME simulation is due to the convention used here for t tr . 
Even though the differences in experimental PR-TRMC 
mobilities are relatively small, we are able to reproduce 
the experimental results not only qualitatively but also 
quantitatively. This suggests that our simulation pro- 
cedure and the underlying morphologies are appropriate 
and that the calculated transport parameters are all in 
the correct order of magnitude. Of particular interest is 
the strong sensitivity of /x to the degree of order in the 
columns, which is itself controlled by the side chain type. 
The highest /i is obtained, in experiment and simulation, 
for the linear side-chain derivatives which also present the 
highest degree of intracolumnar order while the lowest /i 
is obtained for the branched side-chain derivative with 
the least degree of order. The good correspondence be- 
tween simulated and PR-TRMC mobilities arises partly 
from the fact that both simulation and experiment probe 
transport in relatively well aligned columns. To simulate 
ToF experiments, large systems must be studied which 
will require a multiscale ansatz [24| . 

In summary, we have determined charge mobilities of 
several derivatives of the discotic liquid crystal hexaben- 
zocoronene. For the first time, three methods were com- 
bined into one scheme: (i) quantum chemical methods for 
the calculation of molecular electronic structures and re- 
organization energies (ii) molecular dynamics for simula- 
tion of the relative positions and orientations of molecules 
in a columnar mesophase, and (iii) kinetic Monte Carlo 
simulations and Master Equation approach to simulate 
charge transport. Applying this scheme to differently 
substituted HBC derivatives we are able to reproduce 
the trends and magnitudes of mobilities as measured by 
PR-TRMC. This study shows that it is possible to under- 



4 




0.0001 



(b) 



1e-10 1e-09 1e-08 1e-07 1e-06 
time (sec) 



150 



| 100 

o 
o 





(c) 



-6 -4 -2 

log 10 (//eV) 2 



1.5 



o 

C/) 



o 



I" 0.5 h* 
o 



~i 1 r 



~i 1 r 



X 


PR-TRMC 


O 


TOF 


□ 


ME 



J I l_ 



-4- 



j A L 



C 10 C 12 C 14 C 16 C 10-6 PhC 12 



FIG. 3: (a) Simulated ToF hole photo-current transients for 
C12, PI1C12 and Cm-e at T — 300 K and F = 10 5 Vcm" 1 . 
Traces are averaged over 200 snapshots with columns of length 
10. (b) Frequency plots of the logarithm of the transfer inte- 
gral squared. In each case, results for C10-C16 are practically 
indistinguishable from C12 on the scales used, (c) Sum of 
electron and hole mobilities calculated by ToF and by ME, in 
comparison with experimental values. 



stand and reproduce experimental charge transport pa- 
rameters, and, in future, accurately predict them. This 
scheme is the outset for further studies, which will include 
temperature dependencies as well as solid substrates to 
better rationalize the self-organization of these systems. 
By employing multiscale schemes realistic morphologies 
on a semi-macroscopic scale should also become feasible. 
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